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Abstract. The HEGRA system of imaging air Cherenkov telescopes has been used to 
observe the BL Lac object H1426+428 (z = 0.129) for 217.5 hours in 2002. In this data 
set alone, the source is detected at a confidence level of 5.3 cr, confirming this object as a 
TeV source. The overall flux level during the observations in 2002 is found to be a factor 
of « 2.5 lower than during the previous observations by HEGRA in 1999&2000. A new 
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spectral analysis has been carried out, improving the signal-to-noise ratio at the expense 
of a slightly increased systematic uncertainty and reducing the relative energy resolution 
to AE/E < 12 % over a wide range of energies. The new method has also been applied to 
the previously published data set taken in 1999 and 2000, confirming the earlier claim of a 
flattening of the energy spectrum between 1 and 5 TeV. The data set taken in 2002 shows 
again a signal at energies above 1 TeV. We combine the energy spectra as determined by the 
CAT and VERITAS groups with our reanalyzed result of the 1999&2000 data set and apply 
a correction to account for effects of absorption of high energy photons on extragalactic 
background light in the optical to mid infrared band. The shape of the inferred source spec- 
trum is mostly sensitive to the characteristics of the extragalactic background light between 
wavelengths of 1 and 15 yum. 



Key words. BL Lacertae objects: individual H1426+428 - Gamma-rays: observations - 
diffuse radiation - intergalactic medium 



1. Introduction 



The extragalactic objects H1426+428 (Remillard et al 



1989) and Mk n 501 are the only estab 



lished members of the class of extreme synchrotron BL Lac objects dCostamante et al 



2001) 



detected at TeV energies. Generally, these objects are characterized by a flat X-ray spectrum 
with f Y oc v~" with a < 1 extending into the hard X-ray energy band. In the framework of syn- 
chrotron self-Compton (SSC) models, these objects are of particular interest for observations at 
TeV energies, because the inverse Compton peak of the broad band spectral energy distribution 
(SED) is expected to be located at these high energies. The measurement of the position of the 
inverse Compton peak is hampered by the effect of pair-production of high ene rgy photons with 
low energy (optical and infrared) 



Gould & Schreder 



1967 



photons of the extragalactic background light ( Nikishov 



Stecker 



1962; 



de Jagerl) . For the given red-shift z = 0.129 of H1426+428 the 
optical depth exceeds unity even for energies of a few 100 GeV Any detection of a signal at TeV 
energies translates directly into a high luminosity of the source. 

The detec tion and spec t ral measurements of H1426+428 by the CAT, and VERITAS groups 



(Diannati-Atai et al 



2002; 



Petr v et al 



2002) indicate a steep spectrum between 250 GeV and 



1 TeV (well described by a power law dN/dE oc E~" with a photon index of n — 3.66+0.41 for the 
CAT and n = 3.54+0.27 for the VERITAS measurements) whereas at higher energies, a detection 
has been claimed only by the HEGRA group b ased upon observations carried out in 1999 and 



2000 with n - 2.6 ± 0.6 ( Aharonian et al. 



2002a, in the following Paper I) . The spectrum shows 



a flattening at energies above 1 TeV that is consistent with the expected signature of absorption 
as discussed in Paper I. A deep follow-up observation (217.5 hrs) has been carried out in 2002 to 
investigate the energy spectrum with improved statistics. Results of a reanalysis of the original 
data set and the results from the recent observations of 2002 are reported. A possible excess from 
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H1426+428 with a significance of 4.1 <x has been seen with the AIROBICC-Scintillator arrays of 
HEGRA in the data from 1994/95 corresponding to an integral flux of (5.7 + 1.3) • 10 _13 cm _2 s _1 



above E thresh = 21 TeV (Prah 



i: 



1999). The interpretation and the relation of this excess to data 



presented here is not clear at present. 



2. Analysis technique 



The spectral analysis app lied to data fro m the Crab- Nebula ( Aharonia n et al. 



( Aharonian et al 



2 002bL and Mkn 501 (Aharonian el al 



2000), Mkn 421 



1999) benefits from a signal-to-noise 



ratio (S /N) which is usually S /N » 1. Therefore, relaxed cuts (compared to signal search 
optimized cuts) on the mean scaled width (mscw) parameter and the angular separation from 
the source (ff) assures good photon statistics with little systematic effects. However, for weaker 
sources, the signal-to-noise ratio S /N needs to be improved. The relaxed spectral cuts result in 
a poor S IN = 0.1 for H1426+428. A tighter cut on mscw and a smaller angular search bin in- 
creases the S IN to 0.3 reducing at the same time the photon efficiency to the level of 50 %. See 
Table^for a summary of the cuts applied to select the data set, events, and images used in this 
analysis. 

The algorithm for reconstruct ing the energy of the p rimary photon has been improved fol- 



lowing the methods suggested in 



Hofmann et al 



2000). The relative energy resolution for in- 



dividual events has been improved to AE/E < 12 % for a wide range of energies. For energies 
above 1 TeV, the relative energy resolution is as low as 9 %, slightly increasing for lower ener- 
gies. For the threshold region above 500 GeV, the bias of the energy reconstruction method has 
been reduced substantially with respect to the conventional methods. The main contribution for 
this improvement stems from a correction for the position of the shower maximum. Low energy 
air showers triggering the telescopes tend to penetrate deeper into the atmosphere and therefore 
are reconstructed at higher energies. This strong bias is overcome and allows for a more reliable 
energy reconstruction at threshold energies. The r econstruction of the cor e position has been im- 



1999). The improved energy 



proved by assuming a point-like source position ( Hofmann et al. 
resolution reduces spill-over effects distorting the reconstructed energy spectrum in the presence 
of sharp cut-off features in the initial spectrum. 

The scaling of individual image widths to the expectation for simulated showers has been 
modified in order to benefit from the improved core reconstruction. This results in a slightly 
narrower distribution of the mscw-parameter with respect to the conventional core reconstruc- 
tion method. We have cross-checked the background subtracted distribution of mscw for strong 
sources with Monte Carlo predictions and see a very good agreement of the distribution between 
data and simulation. 

The complete analysis chain has been tested with data taken on the Crab-Nebula. The results 
on flux and spectral shape as published in 



Aharonian et al 



(2000) are reproduced. The systematic 



error on the reconstructed photon index is estimated to increase from 0.05 in the case of relaxed 
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Table 1. Listed are data, image- and event-selection criteria. n t ei and nj mage indicating the number 
of active telescopes participating in the data taking and the number of images used for the event 
reconstruction respectively. R meas . and Rqr are the average measured and expected trigger rate 
of the telescopes. r core is the distance of the shower impact point to the position of the central 
telescope. The altitude angle above horizon is indicated as alt, whereas 6 2 is the squared angular 
distance of the shower direction to the source direction. The y-hadron separating quantity mscw 
is further explained in the text. The image parameters distance and size describe the position of 
the shower image centroid with respect to the center of the camera and the total amount of light 
(in photo electrons, p.e.) measured in the image. 



Data selection 



Event selection Image selection 



«td > 3 "image > 2 

1 1 " Rm^./(R C R(0, tt,el)>l W < 200 m 

< 0.2 alt > 45° 

6 2 < 0.014 deg 2 
mscw < 1.1 



distance < 1.7° 
size > 40 p.e. 



Table 2. Listed are for the two data sets: the individual event statistics (N on , N s, and N s = 
N OD - Qnn/Qnff ; Miff), the rate o f excess events per hour (N s ), and the significance S calculated 



using Eq. 17 of Li & Ma 



( 1983). The quantity Q on /Q ff is the normalization for the background 



estimate derived by taking the ratio of the solid angle subtended by the two separate regions. 





reanalyzed 
1999-2000 


2002 


r obs [hrs] 


42.6 


217.5 


(alt) n 


70.8 


68.4 


median(a/f) [°] 


73 


70.4 




309 


1095 


N g 


1447 


6425 


n on /n ff 


1/7 


1/7 




102.3±16.6 


177.1±33.5 


iV s [l/hr] 


2.4±0.4 


0.81 ±0.15 


S[<r] 


6.1 


5.3 



cuts analysis to 0.08 for the analysis presented here. However, given the dominating statistical 
error on low-flux sources, this can be tolerated. 



3. Results 

The data taken on H1426+428 in 1999 and 2000 with the HEGRA system of imaging air 
Cherenkov telescopes have been reanalyzed applying the modified and improved methods de- 
scribed above. The data set has been described in our earlier publication (see Table 1 in Paper I 
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Table 3. Event statistics for individual energy bins. 
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1999-2000 2002 



E 


N 0B 


N oB 


S 


Flux 


N m 


AW 


S 


Flux 


[TeV] 






[o-] 


[10- 13 ph cm- 2 s" 1 TeV" 1 ] 






[0-] 


[10- |3 phcm- 2 s-' TeV" 1 ] 


0.78 


122 


576 


3.78 


187.3 +78.7 


366 


2187 


2.75 


12.3 + 37.8 


1.48 


86 


456 


2.29 


8.11 ±3.51 


394 


2449 


2.16 


2.02+1.40 


2.82 


42 


179 


2.75 


1.92 ± 0.75 


159 


826 


3.33 


0.66 ± 0.24 


5.37 


16 


46 


2.83 


0.47 ± 0.20 


60 


226 


4.00 


0.18 + 0.05 


10.2 


4 


12 


1.36 


0.036 ± 0.034 


7 


58 


-0.43 


< 0.03 (99%c.l.) 



for details on the dates of observations). The signal has been recovered with a slightly higher 
significance than previously claimed increasing the significance from S (Paper I) = 5.8 cr to 
S = 6.1 cr. Note, we have not performed a sensitivity-optimized source search. The cuts are cho- 
sen to perform a spectral analysis (which includes a constraint on the core distance). The deep 
exposure of 217.5 hrs in 2002 confirms the object as a y-ray emitter with a significance for the 
signal of S — 5.3 cr (see Tabled The combined data set with an observation time of 260.1 hrs 
results in S = 7.5 cr (see Fig^for the distribution of the arrival directions). 

The reconstructed energy spectrum of the reanalyzed data set from 1999/2000 agrees well 
with the previously published result (see Fig. Given the quite different approach in event 
selection and reconstruction technique applied in the present analysis, this result demonstrates 
the negligible influence of systematic effects caused by different analysis techniques. 

The recent deep exposure obtained in 2002 shows a decrease in the time-averaged flux by a 
factor of 2.5 with comparison to the observations in 1999 and 2000 (see Fig. [3. The decrease 
in flux observed at high energies is similar to the decrease of the time-averaged X-ray flux as 
measured by the all-sky-monitor (ASM) on-board the Rossi X-ray timing explorer (RXTE). The 
average of the daily rates of the ASM provided by the ASM/RXTE team for the nights of HEGRA 
observations in 1999 and 2000 is 0.33 + 0.08 counts/sec. For the nights of HEGRA observations 
in 2002, the X-ray rate as measured by the ASM drops by a factor of 1 .7 to 0. 19 +0.04 counts/sec. 

4. Discussion 

In order to correct for the effect of absorption of the energetic photons by interactions with the 
extragalactic background light (EBL) in the near (1-3.5 /mi, NIR) to mid (4.5-20 fim, MIR) 
infrared, three largely different models for the spectral energy distribution (SED) of the EBL 
have been used to calculate the opacity t(E, z). 

The calculations of the optical depth r performed here apply a cosmological model with a flat 
cold dark matter (p m = 0.3 p c , p c = 3H^/(^nG) is the critical density) dominated universe with a 
non-vanishing cosmological constant A (pa = 0.7 p c ) and Ho = 60 km s _1 Mpc . Evolutionary 
effects of the EBL have been neglected except for the expansion of the universe. 
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Fig. 1. The number of events reconstructed as a function of the angular separation to the source 
position: After applying a cut on mscw < 1.1 and 9 2 < 0.014 deg 2 (dashed line) to reject back- 
ground events, an excess with a significance of S - 7.5 cr is observed from the direction of 
H1426+428 for the combined data set encompassing 260.1 hrs of observations. The dashed his- 
togram indicates the expected shape of a signal from a point-like source determined from a 
scaled-down signal of the Crab-Nebula. 



The three models are shown in Fig. [3] together with a selection of measurements and lim- 
its on the EBL from the ultra violet to the far infrared wavelength bands. The contribution of 
the far infrared part of the EBL (A > 100 fim) and the cosmic microwave background to the 
op acity is ne gligible, but ha s been included in the calc ulations of r(z, E). The model presented 



by 



Malkan & Stecker (2001); 



de Jaeer & Steckerl (12002) (in the following MS01) is based upon 



a backward-evolution approach with a fast evolving luminosity function between red-shifts 
and 1.3. For the sake of clarity only one of the two alternative models described in MS01 (fast- 
evolution model) is used in the following discussion. Both models follow the same shape in 
the NIR and differ mainly in the absolute level of the radiation (roughly 20 % lower flux in the 
so-called base-line model with respect to t h e fast- evolution model). 



An alternative model ( Prim ack et al 



2001 



in the following PrimOl) invokes a semi- 



analytical approach where the structure formation in a hierarchical scenario is combined with 
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Fig. 2. The differential energy spectra for the different data sets: The reanalyzed data set taken 
in the years 1999 and 2000 reproduces the initially published spectrum with smaller error bars. 
The solid and dashed curves indicating a possible intrinsic source spectrum and the effect of 
absorption are calculated as those shown later in Fig.0J. The data set taken in 2002 is compatible 
with the 1999/2000 data with a flux level reduced by a factor of 2.5 (dot-dashed curve). See also 
Table |3 for the event statistics and differential flux values for the individual energy bins. The 
displayed points are slightly shifted in energy with respect to each other to allow for easier 
reading. 



an initial mass function following a Kennicut-type distribution in order to calculate the SED 
of individual galaxies including realistic dust re-emission models. Following the approach used 
in Paper I, a simple interpolation of the EBL data (labelled Aha02) has been used as a third 
alternative. 

The measure ments and limits on the EBL indicated in Fig. [3] have been taken from 



Hauser & D wek (2001) and updated for recent measurements in the NIR ( Matsumo to et al 
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Fig. 3. The spectral energy distribution (SED) of the of the extragalactic background light (EBL). 
The c ompilation of mea surements and limits between ultraviolet and far infrared are taken from 



Hauser & Dwek 



(2001). Recently publis hed results ha ve been added: the satellite based mea- 



surements between 1 and 4 //m ( Matsumoto et al 



2000), the combined DIRBE/2MASS analyses 



claiming a detection at 1.25 and 2.2 //m (Cambresv, Reach, Beichman, & 



analyses of the 15 /im data taken with the ISO satellite (Elbaz et al 



Jarrett, 



2001), a new 



2002), a model-dependent 



meas urement, and a firm upper-limit at 60 um ba sed upon fluctuation measurements in IRAS 



data (Miville-Deschenes, Lagache, & Puget 



2002). Additionally, three largely di fferent m odels 



for the spectral energy distribution (SE D) of th e EBL are included based upon 



2001) (PrimOl, solid curve), 



Primack et al 



Malkan & Stecker (2001) (MS01 "fast evolution", dashed curve), 



and a phenomenological description of the data (dash-dotted curve, labelled Aha02), interpolat- 
ing between the measured points (Paper I). 



2000 



Cambresv. Reach. Beichman. & 



Miville-Deschenes, Lagache, & Puget 



Jarrettt 120011) . M 
I2OO2L 100, 140, 



2001), MIR (Elbaz et al 



2002), and at 60 Lira 



and 240 fim (Lagac he. Dole. & Pugell 



20031) . Whereas at 60 um the experimental situation is unclear by a factor of » 3 between the dif 



ferent detections ( Finkbeiner . Davis. & Sch lege] 



2000; Miville-Deschenes. Lagache. & Puget 



2002), the claimed values for A = 100 fim are well consistent with each other. 



The opacities Tprfmoi (z, E), TAhatHfo E) have been calculated by integratin g over the respective 
SED w hereas for the MS01 model the prescription of the authors as given in 



de Jager & Stecker 



2002) has been used. 
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Fig. 4. The three panels (a-c) display the observed differential flux values multiplied by E 2 . The 
filled symbols indicate the combined observed spectra, whereas the open symbols indicate the 
intrinsic spectrum inferred from the observed data by applying a model-dependent absorption 
correction. The solid lines are fits of a power law (a) and broken power law (b-c) functions 
to the intrinsic spectrum. The lower dashed curves are the absorption corrected fit functions, 
a) Based upon the PrimOl EBL (see also Fig. O, the absorption corrected spectrum follows a 
power law E 2 dN/dE ocE 05 . b) The fast increase of absorption predicted above 200 GeV for the 
Aha02 EBL leads to a steeply rising (E 2 dN/dE ocE 3 ) source spectrum up to « 700 GeV. Data at 
higher energies are consistent (as claimed in Paper I) with a E 2 dN/dE oc const., c) Applying an 
absorption correction according to the MS01 model: We note a marked upturn of the corrected 
spectrum at energies above « 1 TeV with E 2 dN/dE oc E 3 . Data below energies of 1 TeV follow 
E 2 dN/dE oc const.. 



The wide energy range from 700 GeV up to 10 TeV covered by the HEGRA energy spec- 
trum of H1426+428 offers the unique opportunity to probe the energy region where the source 
spectrum is heavily absorbed by pair-production processes on the NIR/MIR part of the EBL. 

The lower energy spectra published by the CAT and VERITAS collaboration 



(Djannati-Atai et al 




2002) above an energy threshold of « 250 GeV pro- 
vide additional information on the spectral shape at energies where the absorption effect is of 
less importance than for energies above 700 GeV that are accessible by the HEGRA instrument. 
The CAT group has performed a spectral analysis applying a forward-folding technique. The 
residuals of the reconstructed energy spectrum with respect to the forward-folded power law 
assumption contain information on possible deviations from a pure power law. The actual data 
points shown in Fig.|4]are derived by taking the residuals and converting them into differential 
flux measurements. The caveats of using individual data points are clearly larger errors on the 
data points and negligence of the correlation of the energy bins which in turn depends on the 
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assumed spectral shape. However, this method allows to compare the CAT measurements with 
the results of other spectral analysis methods, which suffer from the same caveats. 

The combined data sets for H1426+428 based upon measurements with the CAT, Whipple, 
and HEGRA telescopes (reanalyzed 1999-2000 data set) are presented in Fig.0]a-c. The 2002 
data set has not been considered here, because of the lower flux level during that observation. 
The agreement of the different measurements is very good despite the fact that the observations 
were not carried out simultaneously. The X-ray flux as measured by the ASM averaged over 
the individual observational windows for CAT (0.2 + 0.4 cts/sec), Whipple (0.26 + 0.07 cts/sec), 
and HEGRA (0.33 + 0.08 cts/sec) are in reasonable agreement with each other. The overlap of 
the HEGRA measurements with the results from Whipple and CAT between 700 GeV to 1 TeV 
underlines the consistency of the different measurements. The CAT and VERITAS spectra are 
consistent with each other and no ad hoc normalization is required to combine the different data 
sets. 

After applying a correction to the observed spectrum by multiplying the individual data points 
by exp(Tpri m oi (z = 0.129, E)) we obtain a source spectrum consistent with a hard power law with 
a photon index of 1.5 (see solid and dashed lines in Fig.|^-c respectively for the source spec- 
trum and the observable spectrum). A ^-minimization keeping the photon index fixed at 1.5 
results in a x 1 jd.o.f. = 3/11. This intrinsic source spectrum would imply a position of the in- 
verse Compton peak in the broad-band SED beyond 10 TeV. The energy flux of the unabsorbed 
inverse Compton component exceeds 10~ 10 erg crrT 2 s . An initial modeling with leptonic emis- 
sion models of the broad band SED and specifically the har d source spectrum at TeV energies is 



possible, but requires an additional external radiation field (Costamante et al 



2003). Given the 



low X-ray energy flux observed from this source of a few 10 11 erg cm 2 s 1 during contempo- 
raneous observations, the TeV observations strongly suggest that the y-ray flux dominat es over 



1998) where 



the X-ray flux. This would violate the proclaimed BL Lac sequence ( Fossati et al. 
the energy flux from the synchrotron radiation is expected to be comparable or to dominate over 
the inverse Compton energy flux for TeV blazars. 

The possibly large photon density of the EBL between 1-2 <im as suggested in Paper I to 



accommodate the NIR data by 



Matsumoto et al. 



(2000); 



Cambresv, Reach, Beichman, & Jarrett 



(2001) results in a strong absorption for y energies below 1 TeV. As a consequence of this, the 
inferred source spectrum shown in Fig.|3) rises sharply between 0.25 and 0.7 TeV (E 2 dN/dE oc 
E 3 ) with a subsequent flattening to a constant value. The inferred source spectrum obtained by 
taking HEGRA data alone confirms the result as given in Paper I where the spectral shape of the 
TeV- and X-ray data was tentatively interpreted as self synchrotron-Compton emission. However, 
the inferred rising lower energy part of the unabsorbed source spectrum is not readily explained 
by such a model. 

In the framework of the MS01 model, the flat slope of the SED of the EBL between the NIR 
and the MIR produces a rapidly increasing optical depth for energies above 1 TeV, exceeding 
t = 5 at 2 TeV. Essentially the same behaviour is seen for the somewhat lower EBL level given 
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by the base-line evolution model. After applying the correction for absorption based upon the 
MS01 model, the implied source spectrum is consistent with a power law type source spectrum 
with a photon-index of 2 up to 1 TeV (corresponds to E 2 dN/dE oc const.). Trying to fit the data 
points beyond 1 TeV with E 2 dN/dE oc const, results in a poor x 2 /d.o.f. = 10/3. A broken- 
power law with E 2 dN/dE oc E 3 above 1.4 TeV is a good description of the intrinsic spectrum 
(X 2 /d.o.f. = 3.2/10. A similar upturn is observed for the base-line evolution model (omitted in 
Fig.|4]for the sake of clarity). 

This implied steep upturn in the corrected source spectrum (E 2 dN/dE oc E 1 ) is not easily 
explained by commonly used models invoking either a leptonic or hadronic origin of y-ray emis- 
sion of BL Lac objects. In the case of y-ray production from 7r°-decay from monoenergetic n° a 
E 2 dN/dE « E 2 type spectrum could be produced. However, even in this extreme case, a steep 
rise as it is observed here is not explainable (see e.g. Pohl & Schlickeiser (2000)). 

A sharper pile-up could be expected by speculating that the radiation is the result of bulk 
motion Comptonization of ambien t low-energy thermal photons by a cold conica l wind with 



bulk Doppler factor exceeding 10 7 llAharonian. T imokhin. & Plvasheshnikov 



2002). 



In summary, independent of the details of the model for the EBL, the HEGRA observations 
of H1426+428 between 700 GeV and 10 TeV are strongly affected by absorption of the EBL 
between 1-15 fim. The possible strong absorption combined with the observable flux at TeV 
energies implies an intrinsic spectrum which dominates the broad-band SED of H1426+428. 
This is not in agreement with the conceptual view of the TeV-Blazars being the end-point of 
a phenomenological sequence where the synchrotron component is comparable or dominates 
over the Compton component. The position of the putative Compton peak is possibly at energies 
beyond 10 TeV. As has been shown here, the intrinsic source spectrum strongly depends upon 
the shape of the SED of the EBL. The different models of the EBL may be excluded by requiring 
the intrinsic spectrum to be explainable in the framework of conventional source models. 
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